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The study of b-hadrons made possible a great number of benchmark results in flavour physics. In recent years a lot of 
effort in understanding their dynamics has been done. These studies can represent a possible avenue for the discovery 
of physics beyond the Standard Model (SM) in less studied sectors, and can as well be used as a tool to properly test 
the hadronic calculations reliability. In this sense, a unique opportunity is represented by the B s mesons. They are 
less studied and known with respect to B + and B° mesons, but high precision measurement are possible with the 
current available statistics at b-factories and at hadronic machines. Here we present a brief collection of recent B s 
mesons results. 
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Non-leptonic two-body charmless decays of neutral B mesons (B? s j — > h + h~, where ft, is a charged pion or kaon) 
^ 1 are very interesting for the understanding of flavor physics and CP violation mechanism in the B-meson sector. Their 
^ ■ rich phenomenology offered several opportunities to measure and constrain the parameters of the quark-mixing matrix 
O-i (i. e., Cabibbo-Kobayashi-Maskawa, CKM). These processes also allow to access the phase of the V u b element of the 
CKM matrix (7 angle), and to test the reliability of the SM and hadronic calculations. The presence of New Physics 
(NP) can be revealed by its impact on their decay amplitudes, where new particles may enter in penguin diagrams. 
7—1 ■ The large production cross section of b hadrons of all kinds at the TeVatron and LHC allows extending our 
knowledge of the B® and A° decays, which are important to supplement our understanding of B° and B + meson 
decays provided by the B- factories. The branching fraction of B® — > K~tt + decay mode provides information on the 
CKM angle 7 [1] and the measurement of direct CP asymmetry could be a powerful mo del- independent test of the 
source of CP asymmetry in the B system [2]. The B® — > n + ir~ and B° — > K + K~ decay modes proceed through 
O^n annihilation and exchange topologies, which are currently poorly known and a source of significant uncertainty in 
; many theoretical calculations [3, 4]. A measurement of both decay modes would allow a determination of the strength 
7— I of these amplitudes [5]. Indeed, the CDF Collaborations reported the first evidence for B® — > tt + tt^ decay mode and 
set a limit on the Branching Ratio B(B° — > K + K~) [6]. The LHCb Collaboration confirmed the evidence and also 
J^j " obtained the observation for the B® —> tt + tt~ and set a limit on B(B° —> K + K~) [7]. 

The measurements performed at the two hadron- machines face the same challenge, i. e., to disentangle different 
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decay modes overlapping into a single peak. In this conditions, every mode is a background for the other signals. 
The analysis strategies are similar: the idea is to disentangle the singular components using kinematic and Particle 
Identification (PID) information. While at CDF the PID information allows a limited statistical separation and it 
is necessary combining kinematics and PID information using a 5-dimensional Maximum Likelihood Fit to obtain 
the results, LHCb benefits from the different detector structure and skills. In particular, the presence of the RICH 
detectors allows a powerful particle identification. The efficiency in identifying the final states particles is good 
enough to make possible measurements competitive (and better) with the CDF ones performing only a fit on the 
mass distribution. 



1.1. CDF analysys 

CDF analyzed an integrated luminosity J Cdt ~ 6 fb _1 sample of pairs of oppositely-charged particles with 
Pt > 2 GeV/c and +Pt(2) > 5.5 GeV/c, used to form B candidates. The trigger required also a transverse 

opening angle 20° < A</> < 135° between the two tracks, to reject background from particle pairs within the same 
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jet and from back-to-back jets. In addition, both charged particles were required to originate from a displaced 
vertex with a large impact parameter (100 /im < d ( 1 , 2 ) < 1 mm), while the 6-hadron candidate was required to 
be produced in the primary pp interaction (do < 140 /im) and to have travelled a transverse distance Lt > 200 /im. 
A sample of about 3 million B — » hh! decay modes (where B = B°,B® or A° and h = K or 7r) was reconstructed 
after the off-line confirmation of trigger requirements. In the offline analysis, an unbiased optimization procedure 
determined a tightened selection on track-pairs fit to a common decay vertex. The offline selection is based on a 
more accurate determination of the same quantities used in the trigger. 

No more than one B candidate per event is found after this selection, and a mass (m^) is assigned to each, using 
a charged pion mass assignment for both decay products. 

The resulting 7T7r-mass distributions (see Figure 1, (a)) show a clean signal of B — » hh' decays. Backgrounds include 
mis-reconstructed multibody 6-hadron decays (physics background, causing the enhancement at m WT < 5.16 GeV/c 2 ) 
and random pairs of charged particles (combinatorial background). In spite of a good mass resolution (ss 22 MeV/c 2 ), 
the various B — > hh' modes overlap into an unresolved mass peak near the nominal B° mass, with a width of about 
w 35MeV/c 2 . 



CDF Run II Preliminary J/_ df = 6.11 ftr 1 

□ B°-> K-7C- H KV 
| K*K" 

□ B°— > 71*71- 

"I B°^ K*K" 



K-> P K " 
B?-» rrtt 




Multibody B decays 

I Combinatorial bkg 



5.6 5. 

rrw [GeV/c*] 



(b) 



> 
CD 
in 
CO 



CO 
> 



CO 



32 r 



30- 



« 28 



26- 



CDF Monte Carlo 



A B - 
■ B°- 

□B°- 
o B°- 



■ 7C 7C 

■ KV 

7C*K 

• K*K 




8 -1.0 -0.5 0.0 0.5 1.0 

p = (p -p.)/(p +p.) 



Figure 1: Invariant mass distribution of B — > hh' candidates passing cuts selection, using a pion mass assumption for both 
decay products. Cumulative projections of the Likelihood fit for each mode are overlaid (a). Profile plots of the square 
invariant 7T7r-mass as a function of charged momentum asymmetry /3 for all B° simulated signal modes (b). 



The resolution in invariant mass and in particle identification, provided by specific ionization energy loss (dE/dx) 
in the drift chamber, is not sufficient for separating the individual B — > hh' decay modes on an event-by-event basis, 
therefore a Maximum Likelihood fit, incorporating kinematics and PID information, was performed. The kinematic 
information is summarized by three loosely correlated observables: (a) the square of the invariant 7T7t mass m^; (b) 
the signed momentum imbalance /3 = (p+ — p~)/(p+ +P— ), where p+ (p~) is the the momentum of the positive 
(negative) particle; (c) the scalar sum of particle momenta p to t = P+ + P-- 

The likelihood exploits the kinematic differences among modes (see Figure 1, (b)) by using the correlation between 

the signed momenta of the tracks and the invariant masses rn^ of a candidate for any mass assignment of the decay 

products (m_|_,m_), using the equation 



m, = m 2 — 2rn 2 + mh + m 2 _ + 




This procedure is useful to obtain statistical separation power between 7T7t and KK (or Kit) final states, and 
therefore is a key tool for the measurement of the observables of interest. Kinematic fit templates are extracted from 
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simulation for signal and physics background, while they are extracted from an independent sample for combinatorial 
background. Mass line-shapes are accurately described according for the non Gaussian resolution tails and for the 
effects of the final state radiation of the soft photons. The dE/dx is calibrated over the tracking volume and time 
using about 3.2 millions of D* + — > D°[— > K~h + ]'k + decays, where the sign of the soft pion tags the D° flavor. A 
1.4cr separation is obtained between kaons and pions with p > 2 GeV/c, becoming 2.0a for the couples KK/irir and 
K + ir~ I K~it + . dE/dx templates (signal and background) are extracted from the D° samples used in calibration. 

The dominant contributions to the systematic uncertainty are the uncertainty on the dE/dx calibration and 
parameterization and the uncertainty on the combinatorial background model. An additional systematic uncertainty, 
of the order of 10% has been assessed because of a fit bias, found in the estimate of the relative fraction of the 
B° — > K + K~ decay mode. Other contributions come from trigger efficiencies, physics background shape and 
kinematics, 6-hadron masses and lifetimes. 

The signal yields are calculated from the signal fractions returned by the likelihood fit. For the first time significant 
signal is seen for B° s — > tt + tt~ , with a significance of 3.7a, while the significance for the B° — > K + K~ decay mode is 
2.0 a. 

Absolute results arc listed in Table I; they are obtained by normalizing the data to the world-average of B(B° — > 
K + ir~) [8]. A 90% of confidence level interval is also quoted for the B° — > K + K~ mode. 

The branching fraction of the B° — > K + K~ mode is in agreement with other existing measurements [8], while it 
is higher than the predictions [9] [10]. 

The branching fraction of the B® — > ir + ir~ mode is in agreement with the theoretical expectations within the 
pQCD approach [11], [12] while is higher than most other predictions [9] [13] [14] [15]. 

1.2. LHCb analysis 

We report preliminary measurements on by LHCb collaboration using an integrated luminosity of 0.32 fb _1 col- 
lected during the first part of 2011 at centre of mass energy of 7 TeV. The B9 ■> — > h + h ~ decays events used arc 
extracted from the triggered data using different offline selections, each one targeted to achieve the best sensitivity 
on the measurements of interest. The key point of the analysis is the PID information: the PID observables allows to 
discriminate between the various decay modes. Hence, in order to determine the amount of cross-feed backgrounds 
for a given channel, the relative efficiencies of the PID selection cuts, employed to identify the specific final state of 
interest, play a key role. Huge statistics sample of D* + — > D a (KTr)n + +c.c, decay modes with similar kinematic 
features of the B^ — > h + h ~, are used to extract the PID efficiency. Since the mass region greater then 5.6 GeV/c 2 
are characterized by the presence of A , — > pn(K) decay modes, the PID information has also been studied for protons 
particles using a sample of A — > pir~ decays. 

The B° — > K + K~ and B° s — > tt + tt~ signal events are extracted with an unbinncd maximum likelihood fit to 
different final state mass spectra of events passing the offline selections. The various final states are separated 
using PID requirements in order to have exclusive categories corresponding to distinct final state hypotheses ( 
K + ir~ , K~ir + , K + K~ or 7r + 7r~). In the fits the amount of background for a given channel, due to the other channels 
where at least one particle has been mis-identified (cross- feed background), has been taken into account. The result 
of the fit is shown on Figure 2. 

In particular, LHCb find a significance of more than 5a for the B® tt + tt~ signal, representing the first observation 
of this decay mode. LHCb results are in agreement with CDF results, with about the same level of precision. 



Mode 


Absolute B(10~ 6 ) Absolute Z3(10" 6 ) 
CDF LHCb 


B° -> K+K~ 


0.23 ± 0.10 ± 0.10 0.13 to °ot ± 0.07 
0.57 ± 0.15 ± 0.10 0.98 to'll ± O- 11 



Table I: Absolute branching fractions results for CDF [6] and LHCb [7]. 
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Figure 2: Invariant mass distributions for the mass hypotheses (a, b) K + K~ and (c, d) -k + tt~ . Plots (b) and (d) are the 
same as (a) and (c) respectively, but magnified to focus on the rare B° — > K + K~ and B® —¥ ir + ir~ signals. The results of the 
unbinned maximum likelihood fits are overlaid. The main components contributing to the fit model are also shown. 



2. B° s -> J/^/ (980) lifetime 

In the standard model, the mass and flavor eigenstates of the B® meson differ. This gives rise to particle-antiparticle 
oscillations, which proceed in the SM through weak interaction processes, and whose phenomenology depends on 
the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. The time evolution of B® mesons is governed by a 
Schrodinger equation which contains two 2x2 matrices, called mass and decay matrix. The off-diagonal elements are 
related to observable quantities, namely the mass difference Am s , the phase <f> s between the off-diagonal elements 
of mass and decay matrices, and the decay width difference Ar s , depending on the decay matrix elements and on 
the phase <f) s . A special feature of the B® system is the large value of AT S , which yields a significant difference in 
the lifetimes of the two mass eigenstates of the B®. Within the standard model <fi s is predicted to be very small, 
which in consequence means that CP and mass-eigenstates coincide. If new physics is present, it could enhance </> s to 
large values, a scenario which is not excluded by current experimental constraints. In such a case the correspondence 
between mass and CP eigenstates does not hold anymore and the measured lifetime will correspond to the weighted 
average of the lifetimes of the two mass eigenstates with weights dependent on the size of the CP violating phase <j> s . 
Thus a measurement of the B® lifetime in a final state which is a CP eigenstate provides, in combination with other 
measurements, valuable information on the decay width difference Ar s and the CP violation in B® mixing. 

Using pp collision data with an integrated luminosity of 3.8 fb^ 1 collected by the CDF II detector at the Tevatron 
it is possible to measure the B® lifetime using B® decays to the CP-odd final state J/^>/q(980) with J/ip — ► 
and / (980) -> tt+tt" [16]. 

To extract the B® lifetime CDF used a maximum Likelihood fit on three variables: the invariant mass, the decay 
time and the decay time uncertainty of each candidate. An accurate study of the templates of the signals and of the 
combinatorial background was performed. Figure 3 reports the projection on the life time variable. The dominant 
contributions to the systematic uncertainty are the uncertainties related to the background models, for the mass 
template and for the decay time template. 
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Figure 3: Decay time distribution with fit projection overlaid. 



The result obtained from the fit is: 

r(B° -> J/V)/ (980)) = ( s<a< -) ± °- 03 (syst.)ps. (2) 

This is the first measurement of the B® lifetime in a decay to a CP eigenstate and corresponds in the standard model 
to the lifetime of the heavy B® eigenstate. The measured value agrees well both with the standard model expectation 
as well as with other experimental determinations. While the precision of the lifetime measurement is still limited 
by statistics, it provides an important cross-check on the result determined in B® — > J/ip(j) decays, which relies on an 
angular separation of two CP eigenstates. Furthermore, the measured lifetime can be used as an external constraint 
in the B® — > J /ipcj) analysis to improve the determination of the CP-violating phase. The lifetime measurement in 
B® — > J/ipfo (980) decays is also the next step towards a tagged time dependent CP violation measurement, which 
can provide an independent constraint on the CP violation in the mixing. 
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A B® meson can oscillate into its antiparticle via second order weak interaction transitions, which make its time 
evolution sensitive to contributions from new physics processes. Such contributions are not well constrained yet and 
might be responsible for the deviation from the standard model reported in [17]. The B® eigenstates with defined 
mass and lifetime, B° sL and B® H are linear combinations of the B® and B s states and, in the standard model, 
correspond in good approximation to the even and odd CP eigenstates, respectively. In the absence of substantial 
CP violation, a sizable decay width difference between the light and heavy mass eigenstates, Ar s = T s l — ^ s h, 
arises from the fact that decays to final states of definite CP are only accessible by one of the mass eigenstates. 
The dominant contribution to Ar s is believed to come from the B® —> Dg + Ds (c.c. modes implicitly included), 
which are predominantly CP even and saturate Ar s under certain theoretical assumptions [18, 19], resulting in the 
relation 

AF 

2B(B° S H- « rs + A ^ /2 (3) 

where T s = (T s l + T s h)/2 [20] . However, three-body modes may provide a significant contribution to Ar s [21]. 
A finite value of Ar s improves the experimental sensitivity to CP violation because it allows one to distinguish the 
two mass eigenstates via their decay time distribution. Furthermore, the B® — > Dg Ds decays could be used in 
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future to measure directly the lifetime of the CP -even eigenstate, which would complement the CP odd eigenstate 
lifetime measurement in B® — > J/xpf° (980) as we mentioned before and provide additional information in the search 
for new physics contributions to CP violation. The B® — > Di*^ + decays have been previously studied by the 
ALEPH, CDF, DO, and Belle collaborations [22-25]. Here we present the latest preliminary result from Belle and 
from CDF. 

3.1. Belle results 

At the Belle, the decays B° s ->• D i * )+ D { * ) ~ are reconstructed in a data sample corresponding to an integrated 
luminosity of 121.4 ftr 1 at the T(5s) resonance. The reconstruction of D + candidates has been made using six 
final states: 4>ir+, K°K+, K*°K+, <pp+, K°K*+, and K*°K*+. Belle combined the D+ candidates with photon 
candidates to reconstruct D* s + — > Dfj decays. The strategy of the analysis involve the use of a two dimensional 
Likelihood fit to separate the signal from background. The variables of the fit are the beam-energy-constrained mass 
Mf, c = earn — Pb > an< ^ tne ener gy difference AE = Eb — ^beam' where pb and Eb arc the reconstructed 
momentum and energy of the B® candidate, and -Ek eam is the beam energy. All signal template are extracted from 
the simulations and calibrated using B° -> 7T+ and B° -> D ( * )+ D~ decays. 

The results [27] arc B° s -> D+ Dj = (0.6 ± 0.1 ± 0.1)%, B° -> D*±D*T = (1.8 ± 0.2 ± 0.4)% and B° s -> D*+D*~ = 
(2.0 ± 0.3 ± 0.5)%, the latest corresponding to the first observation. 

Assuming these decay modes saturate decays to CP even final states, the branching fraction determines the relative 
width difference between the CP odd and even eigenstates. Taking CP violation to be negligibly small, we obtain 

AT S /T S = (9.0 ±0.9 ±2.2)% (4) 

where T s is the mean decay width. This result is in good agreement with the current experimental measurements [8] 
and it is consistent with theory [26] . 

3.2. CDF results 

At CDF, the decays B° s -> d£* )+ .d£* ) ~ with D+ -> R-K+TT+ are reconstructed in a data sample corresponding 
to an integrated luminosity of J Cdt ~ 6.8 fb _1 sample collected by the detector. CDF measured the B s production 
rate times the B® — > branching ratio relative to the normalization mode B° — > DfD~. The strategy 

analysis is simple: the relative branching fractions are determined by a simultaneous fit on the mass variable to 
two signal and two normalization samples (fig. 4). The values returned by the fit must be corrected for the relative 
efficiencies. An accurate determination of the efficiencies is achieved by taking into account the Dalitz structure of 
the D+ decay. As a result CDF obtained [28]: B(B° S -> D+D~) = (0.49 ± 0.06 ± 0.05 ± 0.08)%, B{B° S -> D^D**) 
= (1.13 ± 0.12 ± 0.09 ± 0.19)%, B(B° S -> D*+D*~) = (1.75 ± 0.19 ± 0.27 ± 0.29)% and B(B° S -> L>i* )+ L>i* )_ ) = 
(3.38±0.25±0.30±0.56)%. Statistical, systematic and normalization uncertainties are reported. These measurements 
represents the world best measurements up to date. 

The information on the branching ratios are useful to infer indirect information about AT Sy using relation 3: 

Aryr s = (6.99 ± 0.54 ± o.64 ± 1.20)% (5) 

which is consistent with SM expectation [26]. 

4. CONCLUSIONS 

These measurements represent a brief collection of different attempts to enter a new era: with the current statistics 
from b-factories and hadron machines, it is possible to undergo the systematic and detailed exploration of the B s 
sector. Many new results are expected to come: Belle and CDF are squeezing up their full data sample, while the 
advent of LHCb will make possible unprecedented statistical precision measurements. 
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Figure 4: Invariant mass distribution of (a) B° -> Z)+( (? !>7r + )D s "(< ? !>7r"), (b) S° D+(^*°7<' + )D s "(( ? !>7r"), (c) B° -> 
-D^"(07r + )/? _ (_R' + 7r _ 7r + ), and (d) B° —> Df (K*° K + )D~ (K + ir~ir + ) candidates with the simultaneous fit projection overlaid. 
The broader structures stem from decays where the photon or ir° from the D*^ decay is not reconstructed. Misreconstructed 
signal events in (c) show up as reflections in (b). 
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